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ABSTRACT 

Context. The hot gas in clusters and groups of galaxies is continuously being enriched with metals from supemovae and stars. It is 
well established that the enrichment of the gas with elements from oxygen to iron is mainly caused by supernova explosions. The ori- 
gins of nitrogen and carbon are still being debated. Possible candidates include massive, metal-rich stars, early generations of massive 
stars, intermediate or low mass stars and Asymptotic Giant Branch (AGB) stars. 

Aims. In this paper we accurately determine the metal abundances of the gas in the groups of galaxies NGC 5044 and NGC 5813, and 
discuss the nature of the objects that create these metals. We mainly focus on carbon and nitrogen. 

Methods. We use spatially-resolved high-resolution X-ray spectroscopy from XMM-Newton. For the spectral fitting, multi- 
temperature hot gas models are used. 

Results. The abundance ratios of carbon over oxygen and nitrogen over oxygen that we find are high compared to the ratios in the 
stars in the disk of our Galaxy. The oxygen and nitrogen abundances we derive are similar to what was found in earlier work on 
other giant ellipticals in comparable environments. We show that the iron abundances in both our sources have a gradient along the 
cross-dispersion direction of the Reflection Grating Spectrometer (RGS). 

Conclusions. We conclude that it is unlikely that the creation of nitrogen and carbon takes place in massive stars, which end their lives 
as core-collapse supemovae, enriching the medium with oxygen because oxygen should then also be enhanced. Therefore we favour 
low-and intermediate mass stars as sources of these elements. The abundances in the hot gas can best be explained by a 30-40% 
contribution of type la supemovae based on the measured oxygen and iron abundances and under the assumption of a Salpeter Initial 
Mass Function (IMF). 

Key words. Galaxies: Abundances - X-rays: Galaxies - Galaxies: Groups: Individual: NGC 5044, NGC 5813 - Galaxies: Groups: 
General - Galaxies: Clusters: Intracluster Medium 



. 1. Introduction constrain supernova model s and ratios (see e.g. the review by 

' , . , , ., TT . ... , tn „ i [Bohringer & Werner 2010). While earlier work focussed on 

Most of the galaxies in the Universe reside in groups (e.g. ttullyj du more r£centl ^ metal contents of groups has been 



[1982b. Groups of galaxies typically contain up to tens of galax- investigated in more detail ( Rasmussen & Ponman 2007; 

^- . les. Due to the fact that their gravitational potential is relatively Wemer g ^ 2Q09) Measurements with the RGS are the only 

O ■ shallow, groups of galaxies are good laborat on es tc .study non - tQ constrain ^ carbon and ^ Q abundance in the hot 
, gravitational effects (see e.g. the review by IMulchaevI 120041) . 
I • The gravitational compression and infall shocks, which heat up 

> the gas, depend on the depth of gravitational potential, therefore The abundances observed in the intracluster medium are the 

rj ! groups of galaxies are cooler than the more massive clusters. result of the production of the elements in stars and the trans- 

Since the temperature of groups is rather low (~ 10 7 K), the P ort of the elements released to their current location.The met- 

b ' emission lines of metals present in the medium are relatively als with even atomic numbers from oxygen to up to iron and 

■' strong. This yields an insight in the abundance distribution of nicke l are produced in supemovae: oxygen, neon and magne- 

metals in the intra-group medium. Since all the metals in the sium mainl y in core-collapse supemovae, iron and nickel mainly 

intra-group medium have been produced by stars, metal abun- in tv P e la supemovae. Both types create significant amounts 

dances in the intra-group medium are a tracer of the history of of silicon, sulfur, argon and calcium. The relative fraction of 

stellar evolution. these element s is however unknown and highly mode l depen- 

X-ray spectroscopy of groups and clusters of galaxies is a de nt (see e.g.|lwamoto et al.||l999|;[Woosley et al.||2002|). Carbon 

powerful tool to measure the abundances in the hot gas. The and nitrogen are produced by a variety of sources, the relative 

broad spectral band of the Europ ean Photon Ima ging Cameras importance of which is still under debate (for an overview of 

(EPIC -MOS, iTurner et al.1 l200ll and EPIC-pn. Isiruder et al.1 possible sources of several elements, see |Romano etaLj|2010|, 

2001) and the high s pectral resolution of the Reflection Grating and references therein). Based on optical observations of stars 

Spectrometers (RGS. Iden Herder et al.l200lT) on board of XMM- in ° ur galaxies, most authors conclude that superwinds of mas- 

Newton allow us to determine the key abundances needed to sive metal-poor stars enriched the gas with carbon early in the 

Galaxy disk's evolution. Some authors argue that at a later stage, 

Send offprint requests to: y.g.grange@sron.nl a significant amount of carbon is created by low-mass stars 
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(e.g.lCescutti et alJl2 009). Others fe.g. lBensbv & Feltzingf2 006: 
Mattsson 2010) advocate the enrichment at later stages is mainly 
due to intermediate mass stars, during their As ymptotic Giant 
Branch (AGB) phase. Bensbv & Feltzi ng (2006) conclude from 
the literature that it is likely that the enrichment with carbon is 
probably due to a complicated, and finely tuned set of objects. 
Nitrogen is created by the conversion of carbon and oxygen dur- 
ing the CNO c ycles of hydrogen burning in low and in termediate 
mass stars dEdmunds & Pagell 1978t ICh iappini et al. 2003 ) and 
metal-poor massive stars (Matteucci 1986; Mevnet & Maeder 
l2002h . 

Several mechanisms were proposed that enrich the gas in 
groups and clusters of galaxies with the metals created inside 
the galaxies. An overview of these mechanisms is given by 
Schindler & Diaferio (2008). Simulations suggest that the gas 
is mainly enriched by galactic winds, powered by core-collapse 
supernovae at early stages of the evolution. During the life time 
of the group, the winds get more and more driven by type la su- 
pernovae. Other mechanisms include e.g. uplift of metals by the 
central Active Galactic Nucleus (AGN) and ram-pressure strip- 
ping of infalling galaxies. The relative importance of these pro- 
cesses is however yet unknown. 

In this paper, we study two X-ray bright groups of galaxies, 
NGC 5044 and NGC 5813. Both groups are named after the gi- 
ant elliptical galaxy in their core. Multiple dwarf galaxies are 
found around the central elliptical galaxy of these groups. NGC 
5044 has prev iously been observed with se veral X-ray satel- 
lites (E instein: iFabbiano et al.| [l992l ROSAT: IPavid et al.ll 1994 
ASCA: Fukaza wa et al] 1 19961) . The X-ray emission is sharply 
peaked in the central region and the temperature declines to- 
wards the core. Recently, the kine matic struc t ure of NGC 5044 
was studied with Chandra by IPavid et aD d2009h . Chandra 
images show multiple filame nts and cavities, whic h indicate 
activity of the central AGN dGastaldello et all 120091) . XMM- 
New ton observations sho w a clear multi-temperature structure 
and iTamura et al.l (120031) put constraints on the metal abun- 
dances. Since then, the models and instrumental calibration have 
been improved and longer observations have been performed to 
measure the abundances with a higher precision. 

The group around NGC 5813 is part of a larger gravitation- 
ally bound system which also includes the group around NGC 
5846 . NGC 5813 has a kilop arsec size, kinematically decoupled 
core (Efstathiou et ali i 19821) caused by a histor y of minor or in- 
termed iate mergers long ago. It is classified by Ems ellem et alj 
(12001) as being a slow rotator, which they conclude is an ex- 
treme evolutionary end point reached in a deep potential well. 
In essence, this me ans that this group is extremely relaxed. 
iRandall et al.l d201 ll) studied the shocks caused by the AGN 
feedback in this group using data taken by Chandra, and con- 
clude that the outflow of the central AGN is energetic enough to 
heat up the gas. 

We have obtained deep XMM-Newton observations of these 
two groups. Exploiting the high spectral resolution and sensi- 
tivity of XMM-Newton, we aim to accurately determine metal 
abundances in the hot gas. Especially the carbon and nitro- 
gen abundances contain clues on the likely candidates for their 
production. In addition, the ratio between type la and core- 
collapse supernovae can be constrained using the oxygen and 
iron abundance. We use the measured abundances in the intra- 
group medium to test supernova yield models. 

Throughout this paper, we use Ho = 70 km s _I Mpc -1 , 
Q M = 0.3, and Q A = 0.7. At the redshift of NGC 5044 
(z = 0.009), 1' represents a distance of 11.1 kpc. At the red- 
shift of NGC 5813 (z = 0.007), 1' represents a distance of 



Table 1. Log of the observations used. 





Exposure 


Exposure Exposure 




Position 


obsID 


MOS 


pn 


RGS 


Start date 


angle 




(ks) 


(ks) 


(ks) 


(Y-M-D) 


(degrees)" 


NGC 5044 


0037950101 


21 


12 


17 


2001-01-12 


110 


0554680101 


104 


86 


95 


2008-12-27 


114 


Total exposure 


125 


98 


112 






NGC 5813 


0302460101 


30 


22 


25 


2005-07-23 


291 


0554680201 


60 


44 


51 


2009-02-11 


104 


0554680301 


60 


43 


50 


2009-02-17 


102 


Total exposure 


150 


109 


126 







Notes. Quoted exposures are after filtering the soft-proton background. 
<n) The zeropoint of the position angle is the celestial north. The positive 
position-angle direction is towards the east. 



7.92 kp c. The abundan ces are in units of the proto-Solar abun- 
dances (Lodders 2003). Quoted statistical uncertainties are at a 
l<x (68%) confidence level. Resul ts obtained from the liter ature 
using different abundanc e scales dAnders & Grev esse 1989] and 
Grevesse & Sauval 1998) or significance levels (e.g. 90%) were 
converted to the units of this work before comparing them. 

2. Data reduction and analysis 

All the data used for this research were obtained using the in- 
struments on board of the XMM-Newton satellite. A log of the 
observations is presented in Table Q] The recent (2008 and 2009) 
data were obtained using the thin filter and for the archival data 
the medium filter was used. We use data from both the EPIC 
and the RGS. The EPIC-MOS was in full frame mode during 
all observations. For the EPIC-pn, the extended full frame mode 
was turned on during the 2005 observation of NGC 5813, while 
for all others it was in the Full Frame mode. In total, 153 ks of 
data are available for NGC 5044 and 173 ks for NGC 5813. The 
reduction was done with the XMM-Newton Science Analysis 
Software (SAS) version 9.0.0 with the calibration files created 
on 22 June 2009. 

In order to reduce contamination from soft-proton flares, we 
created a lightcurve for each instrument and observation. For 
EPIC, the light curve was based on the data in the 10-12 keV 
energy band and for RGS we used the data from CCD number 
9, where hardly any emission from the group is expected. We 
binned the light curve in 200 s intervals for RGS and 100 s inter- 
vals for EPIC. A Poissonian distribution was fitted to the count- 
rate histogram. We rejected all time bins for which the number 
of counts lies outside the interval fj. + 2 yfjj, where /i is the fitted 
average of the Poissonian. The clean exposure times are listed in 
TableQ] Since we cut off all data above yu + 2<x as well as all data 
under fi-2cr, we do not expect the flux to change due to our filter- 
ing. To check whether a less severe filtering would significantly 
improve the statistics, we fitted RGS spectra for which we did 
not use any filtering at all. These fits are statistically comparable 
to the results of this work. In addition, we created a light curve 
from the events between 0.3 and 1 keV to check whether our 
spectra are contaminated by charge-exchange emission from the 
interactions between the solar wind and gas in the ecliptic plane 
of our solar system and the Earth's atmosphere. This emission 
should show up as flares in the light curve under 1 keV. The 
light curve does not show significant flares in this energy region. 

The RGS is a slitless instrument. This means that the ze- 
ropoint of the wavelength calibration is dependent on the po- 
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Fig. 1. Smoothed Chandra image of NGC 5044. The edges of 
the RGS for obsID 0037950101 are indicated with black lines. 
The black arrow has length 1', and points in the positive cross- 
dispersion direction. 




Fig. 2. Smoothed Chandra image of NGC 5813. The edges of 
the RGS for obsID 0554680201 are indicated with black lines. 
The black arrow has length 1', and points in the positive cross- 
dispersion direction. 



sition of the source in the detector field of view. We gauge 
our wavelength scale of RGS on the peak of the emission in 
the X-ray image. We also use this position as the source cen- 
tre for the EPIC extraction region. For NGC 5044, we take 
a=13 h 15 m 24 s , £=-16° 23' 08" and for NGC 5813, these coor- 
dinates are a=15 h 01 m 1 I s , (5=01° 42' 06". These coordinates are 
within 2" from the optical position of the central galaxy. 

2.1. Reduction of the RGS data 

The RGS data were reduced using the standard pipeline process- 
ing of the SAS. We subtract from the total spectrum a model 
background spectrum which is created by the standard RGS 
pipeline. This method consists of selecting a template back- 
ground file, based on the count rate in CCD 9 of the RGS. We do 
not explicitly model the cosmic X-ray background in RGS be- 
cause any diffuse line source would be smeared out into a broad, 



nearly flat continuum-like component. We included 100% of the 
cross-dispersion direction in the spectrum because the extent of 
the sources is larger than the width of the RGS. All other options 
were left to their default values. To make one-dimensional spa- 
tial profiles of the parameters that we measured using RGS, we 
also created spectra based on five narrow regions in the cross - 
dispersion direction of the RGS, within 2' from the core. During 
one of the observations of NGC 5813 (obsID 0302460101) the 
position angle differs by ~180° from the other two observations. 
We fitted the spectra from the cross-dispersion strips of this ob- 
servation simultaneously with the spectra from the opposing re- 
gions of the other two observations. In the other cases, we as- 
sumed that the differences in roll angle were not large enough to 
introduce systematic effects. The relative position of the sources 
to the RGS detectors and the positive cross-dispersion directions 
are shown in Figs. Q] and [2] 

Since, as mentioned before, RGS is a slitless instrument, all 
photons within the field of view end up in the final spectrum. 
This causes a broadening of the spectral lines. If a photon origi- 
nates from an angle A6 (in arc minutes) projected along the dis- 
persion direction, it will be registered as having a wavelength 
which is shifted by AA (in A) with respect to its original wave- 
length. The XMM-Newton Users Handbook gives the relation 

A,i = 0.138Al9/m. (1) 

Here, m is the (absolute value of the) spectral order. We ex- 
tracted the brightness profile of the source in the dispersion di- 
rection from the EPIC-MOS image and used this t o account for 
the broadening following the method developed bv lTamura et alJ 
(12004 . We use the spectral data between 7 A and 36 A. This up- 
per limit is chosen so that we include the C vi line at 33.7 A. 

2.2. EPIC data reduction 

Because we want to compare our EPIC results with RGS, we 
derive EPIC spectra from a circle with radius 3 ' around the cen- 
tre. The diameter corresponding to this radius is comparable to 
the cross-dispersion width of the RGS (5')- We identified point 
sources, based on visual selection, and excised a region with a 
radius of 15" around them. In total, we excised one point source 
from the NGC 5044 observation. Two point sources were re- 
moved which lie on the edge of the 3' circle. From the field 
around NGC 5813, three point sources were identified and re- 
moved. Based on the spectrum, we created files containing the 
effective area and response matrix using the default SAS-tasks. 
To correct for vignetting effects, we created these files using a 
detector map of the source. Furthermore, we corrected the pn ob- 
servations for Out Of Time (OOT) events. We ignore all events 
corresponding to energies higher than 10 keV. For EPIC-MOS, 
we ignore the events below 0.4 keV. We also ignore all events 
below 1 .5 keV in the EPIC-pn data because these data are af- 
fected by a gain shift, leading to large-scale systematic residu- 
als below this energy. This effect seems to be present in both 
sources. 

2.2.1. EPIC background 

To correctly fit the E PIC spectra, w e need to subtract or 
model the background. Snowde n et al.l d2004l) used a method 
to subtract the EPIC-MOS instrumental background by scal- 
ing closed filter observations, as created by the XMM- 
Newton Science Operations Centre (SOC). Other authors 
(lArnaud et al.l200lURead & Ponmanl200llDe Luca & Molendil 
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Table 2. Instrumental lines modeled for EPIC-MOS and EPIC- 
pn. 



Element 


Energy 
( keV) 


MOS 


pn 


AlKff 


1.486 


+ 


+ 


AlKjS 


1.557 


+ 


+ 


SiKa 


1.740 


+ 


- 


Si KjS 


1.835 


+ 


- 


TiKa 


4.51 


- 


+ 


CrKa 


5.41 


+ 


+ 


Mn Ka 


5.89 


+ 




FeKa 


6.40 


+ 




NiKa 


7.47 




+ 


CuKar 


8.04 




+ 


CuKjS 


8.90 




+ 


ZnKo- 


8.63 


+ 


+ 


ZnKjS 


9.57 




+ 


Au La 


9.72 


+ 


+ 



Notes. Lines from Ide Plaa et al.l d2010l) and iLec cardi & Molendi 
(2008). A "+" indicates the line is present in the background model 
for the instrument. A "— " means the line is not part of the background 
model for the given instrument. 

2004; iNevalainen et al.l 12005) developed and used comparable 
methods. The methods mentioned are highly dependent on the 
region of the CCD outside the field of view. In the pn, this re- 
gion is very small which makes the method sensitive to sta- 
tistical fluctuations. Since we want to simultaneously fit the 
three instruments, we need to use a unified background subtrac- 
tion method. Our method involves modellin g the background by 
use o f additional spectra l models (see e.g. ILeccardi & MolendU 
2008: Ide Plaa et alJ2010h . as follows. We model both instrumen- 
tal and Cosmic X-ray background. 

For the instrumental background, we used the closed filter 
observations, provided by the XMM-Newton SOC. We extracted 
a spectrum from these closed filter files, using the same proce- 
dure as for the actual data. We modeled this instrumental back- 
ground using a broken power law, which describes the hard par- 
ticle spectrum, and a number of delta functions that represent 
the fluorescence lines caused by the interaction of hard particles 
with the instrument itself. We list the fluorescence lines that we 
take into account in Table [2] Since high-energy particles are not 
focussed by the telescope in the same way as X-ray photons, 
we do not fold these models through the Ancillary Response 
File (ARF), which describes the effective area of the instrument. 
The instrument degrades with time, leading to an increase of the 
instrumental background. In addition, the instrumental power- 
law does not vary over the field of view, while the strength of 
the instrumental lines does depend on the position on the detec- 
tor. During the fits to the core spectra, we fitted the instrumen- 
tal background for each observation by allowing the normalisa- 
tions of the instrumental background components to vary within 
a factor of 4 of the best-fit value, derived from the closed filter 
spectrum. We also left free the power-law break strength. We as- 
sumed that the two observations of NGC 58 1 3 on February 2009 
had the same instrumental background. 

The second group of components models the emission by the 
astrophysical back- and foreground. These X-rays are focussed 
by the mirror in the same way as source photons, and therefore 
these components should be folded through the ARF. At high en- 
ergies (> 2keV), the background is dominated by the Cosmic X- 
ray Background (CXB) of unresolved point sources. This com- 
ponent is modeled using a power-law, the Extra-galactic Power 



Table 3. Relevant parameters of the fitted background compo- 
nents. 



Component 


Parameter 


JNCjC j 044 




PPT 


v a 
i 


1 zLI 


1 zll 




Flux"-* (2-10 keV) 


2.24 


2.24 




Flux a 6 c (0. 1-2.5 keV) 


0.96 


1.02 


LHB 


kT ( keV) 


0.11 


0.07 




Flux* (2-10 keV) 


1.57 x 10~ 7 


2.05 x 10~ 12 




Flux* (0. 1-2.5 keV) 


3.61 


4.48 


halo 


kT ( keV) 


0.29 


0.19 




Flux* (2-10 keV) 


1.84 x 10- 3 


1.23 x 10~ 4 




Flux* (0.1-2.5 keV) 


0.78 


0.90 



Notes. (a) Fixed during fit ib) In units of xl0~ 14 Wnr 2 deg~ 2 ic) Values 
differ because of the difference in galactic H i values. 



Law (EPL). We fix the photon ind ex T to 1 .41 and the 2-10 k eV 
flux to 2.24 x 10~ 14 Wirr 2 deg" 2 dDe Luca & Molendill2004T) . 

At lower energies (< 2 keV), foreground emission from the 
Local Hot Bubble (LHB) dominates the background spectrum. 
Contrary to the other X-ray background components, this one 
is not absorbed by the neutral gas of our Galaxy. Furthermore, 
the diffuse halo or disk of the G alaxy also contributes t o the 
foreground at these low energies dKuntz & Snowdenll2001l) . We 
modeled these two components using a Collisionally Ionisation 
Equilibrium (CIE) plasma model. Since the temperatures of 
these components vary across the sky, we have to fit the tem- 
peratures locally. Therefore, we extracted a spectrum from an 
annulus between 9' and 12' centred on the emission peak of the 
group. In this annulus the group emission is low. In addition to 
the background components, we added an extra CIE component 
to model any residual group emission. The best-fit values of the 
relevant parameters are listed in Table|3] We keep the parameters 
of the X-ray background model components frozen to the best- 
fit values while fitting the spectra in the core. The normalisations 
of the components are scaled to the size of the extraction region. 

2.3. Spectral models 

We fit our spectra using the SPEX package dKaastra et al.ll 996). 
Since clusters and groups are not isothermal, we use multi- 
temperature models. We fit the spectra using two different empir- 
ical parametrisations of a Differential Emission Measure (DEM) 
model. The emission measure Y = n e «hdV with n e and «h the 
electron and proton densities and V the volume of the source. 

The first parametrisation is the gdem model, in which the dif- 
ferential emission m easure dY/dT as a function of temperature 
is distributed as (e.g. Ide Plaa et al.ll2006h 

§ = (2) 
dI o- T V27T 

In this equation, x = log(T) and xq = log(7o). 

The second D EM model we fit is the wdem model 
dKaastra et al.l2 004). The model consists of a number of thermal 
components distributed as a truncated power law. The model is 
characterised by: 

dY f cT l ' a pT m . dx < T < T m , x 

dT \o r>r max vr<ysr max • w 

For a — > 0, the model becomes isothermal with temperature 
kT 

max- 
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Wavelength (A) 

Fig. 3. Fluxed, background-subtracted RGS spectrum of NGC 
5044. The (red) line through the data points is the best fit. The 
most important spectral lines are indicated in the plot. 



CO 




CO I . i .... i .... i .... i 

° -2 -1 1 2 

r (') 

Fig. 4. Iron and temperature profiles in the cross-dispersion di- 
rection based the fits to the RGS data of NGC 5044 using the 
gdem model. 



In order to compare the multi-temperature models, we 
calculate the emi ssion-weighted temperature for wdem (see 
Ide Plaa e~ al. 2006, for the derivation), 



kT n 



(l + l/q)(l-^ +2 ) 1T , 

(2 +l/ ff )(l -/?!/«+!) ' 



(4) 



For gdem the peak temperature is fitted. 

The hot gas models allow us to determine the abundances of 
metals in the intra-group medium. Since the neon and magne- 
sium abundance could not be constrained using EPIC, we fixed 
the ratio between these and iron to the values found in RGS. 
For the Galactic foreground absorption due to neutral gas, we 
used the Hi column density from the Leiden/Argentine/Bonn 
survey dKalberla et all [2005). For NGC 5044 we use a col- 
umn density of 4.87 X 10 24 nr 2 and for NGC 5813 the value is 
4.37 x 10 24 nr 2 . 



3. Results 

We fitted the spectra of NGC 5044 and NGC 5813 using the 
models described in section |2~31 In order to easily compare the 
abundances found to other groups, we quote abundance ratios 
for most elements. We chose to quote the abundances relative to 
iron, because iron has the smallest statistical uncertainty of all 
measured abundances. 



3.1. NGC 5044 

In Fig. [3] we show the RGS spectrum of NGC 5044. The best- 
fit values are listed in Table |4] The spectrum shows a wealth of 
lines that can be used to constrain the abundances and multi- 
temperature structure of the gas. Most lines are fitted well, how- 
ever both the Fe xvm and O vm-LyyS line blend around 16 A and 
the blend of Fe xvm and Fe xix near 14 A are not fitted optimally. 
Near 1 1 A the model is also underestimating the data. 

Table [4] shows that this source has a significant temperature 
structure to which the gdem model fits slightly better than the 
wdem model. The derived abundance ratios, as reported in Table 
[U are consistent. The absolute iron abundances derived from 
these models differ by about 3<x. The mean temperature is con- 
sistent between the two models. 




Fig. 5. O/Fe and N/Fe profiles (in proto-Solar units by lLoddersI 
2003) based on the fits to the RGS data of NGC 5044 using the 
gdem model. The nitrogen abundance could not be constrained 
in the rightmost bin and is therefore not present in the plot. 



The background for RGS at long wavelength is high for 
extended sources, and hence some systematic residuals remain 
near the C vi line at ~ 34 A. 

Fig. [4] shows the profiles of the gdem peak temperature and 
iron abundance as measured along the cross-dispersion direc- 
tion of RGS (see Fig. Q] for the physical direction of the cross- 
dispersion direction). The temperature shows an increase from 
the centre of the galaxy towards the outer parts, as is normal in 
the cooling cores of groups and clusters of galaxies. The iron 
abundance shows a hint of a linear decrease over the cross- 
dispersion direction. Fitting a linear model to the data yields a 
normalisation of 0.663 ± 0.016 solar at r — and the slope of the 
line is -0.08 ± 0.02 solar per arc minute. The profiles of the rel- 
ative abundances of oxygen and nitrogen are presented in Fig. [5] 
The oxygen over iron abundance shows no significant variation 
over the cross-dispersion direction. The nitrogen abundance is 
poorly constrained. The EPIC spectrum, (Fig. |6) shows system- 
atic residuals above 2 keV, where the background is becoming 
more important possibly to a low contribution of quiescent soft- 
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1 2 

Energy (keV) 

Fig. 6. EPIC spectrum of NGC 5044. The top panel shows the 
data points with the model as a histogram through them. Upper 
data set: EPIC-pn, bottom data set: EPIC-MOS. The bottom 
plot shows the fitting residuals. Squares: EPIC-pn, circles:EPIC- 
MOS. 

Table 4. Fitted values to the NGC 5044 spectra. 



Parameter 


RGS wdem 


RGS gdem 


EPIC gdem 


Y(10" nT J ) 


3.02(6) 


2.95(5) 


2.555(10) 


r peak (keV) 


1.020(10) 


0.910(5) 


0.9585(14) 


Tmean (keV) 


0.887(10) 






C7"t 




0.100(3) 


0.1464(12) 


l/a 


0.178(9) 








0.4 






C/Fe 


1.5(3) 


1.5(3) 




N/Fe 


1.4(2) 


1.5(2) 




O/Fe 


0.61(3) 


0.60(3) 


0.472(17) 


Ne/Fe 


0.35(12) 


0.37(10) 




Mg/Fe 


0.80(6) 


0.73(5) 




Si/Fe 






0.808(12) 


S/Fe 






0.871(18) 


Ar/Fe 






0.789(13) 


Ca/Fe 






1.84(13) 


Fe 


0.572(14) 


0.616(16) 


0.782(6) 


Ni/Fe 






1.60(5) 


C-stat / d.o.f. 


2591 / 1332 


2418/ 1332 


3781 / 1775 



Notes . Abundance ratios are given in proto-Solar units by Lodders 
(2003). Numbers between parentheses are the Icr errors on the same 
number of least significant digits. <a) Kept fixed during the fit. 



protons. The residuals also show some systematic effects in the 
iron-L complex at energies below 2 keV. 

The iron abundance derived using the EPIC data is signif- 
icantly higher than the one derived from RGS. For the oxygen 
over iron ratio, the value derived from EPIC is lower than the 
RGS value. The abundance ratios of silicon, sulfur, argon, cal- 
cium and nickel over iron are high. The EPIC fits give a tem- 
perature structure with a higher peak temperature and a lower 
emission measure than for those of the RGS. For the EPIC spec- 
tra, the derived temperature structure is broader than what we 
find for the RGS spectra. 

3.2. NGC 5813 

We present the RGS spectrum of NGC 5813 in Fig.0 The model 
is underestimating the line flux of the oxygen and iron line blend 



E 

<n 
O 




15 20 25 

Wavelength (A) 



Fig. 7. Fluxed RGS spectrum of NGC 5813. The (red) line 
through the data points is the best fit. The most important spec- 
tral lines are indicated in the plot. 




Fig. 8. Iron and temperature profiles in the cross-dispersion di- 
rection based on the fits to the RGS data of NGC 5813 using the 
gdem model. 



near 16 A, which is also the case for NGC 5044. Contrary to the 
fit for NGC 5044, the 14 A blend is well fitted by the model. 

The values and accuracies of the carbon and nitrogen abun- 
dances are well defined and comparable to what we find for NGC 
5044. Also for this source, the background at ~ 34 A introduces 
extra uncertainty in the measured carbon abundance. 

The fits to the RGS spectrum of NGC 5813 show no signifi- 
cant difference between wdem and gdem. Both models converge 
to the single-temperature limit with essentially the same values 
for all parameters and the same statistical significance. 

The RGS profiles along the cross-dispersion direction of the 
gdem peak temperature and the iron abundance that we obtained 
from the RGS data are shown in Fig. [8] For the physical mean- 
ing of the cross-dispersion direction, see Fig. [2] A decrease to- 
wards the centre of the temperature is observed, like in cool -core 
groups. However, the central temperature gradient is shallower 
than for NGC 5044. The iron abundance appears to be asym- 
metrically distributed around the core. A linear fit to this data 
set yields a normalisation of 0.650 + 0.018 solar at r — and a 
slope of -0.13 + 0.03 solar per arc minute. In Fig. |9] the radial 
profiles of oxygen and nitrogen over iron abundances are shown. 
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Table 5. Fitted values to the NGC 5813 spectra. 





r (') 



Fig. 9 . O/Fe 
2003) based 
gdem model 



and N/Fe profiles (in proto-Solar units by lLoddersI 
the fits to the RGS data of NGC 5813 using the 




I 0.5 1 2 5 10 

Energy (keV) 

Fig. 10. EPIC spectrum of NGC 5813. The top panel shows the 
data points with the model as a histogram through them. Upper 
data set: EPIC-pn, bottom data set: EPIC-MOS. The bottom plot 
shows the fitting residuals. Squares: EPIC-pn, circles: EPIC- 
MOS. 



The O/Fe abundance ratio is flat, with a hint of a decrease from 
the negative to the positive cross-dispersion direction. 

The EPIC spectra of NGC 5813 are shown in Fig. [TU] In the 
iron-L complex, around 1.2 keV some residuals remain which 
do not show up as strong in the spectra of NGC 5044. Between 
3 keV and 5 ke V, the systematic residuals are similar to the same 
spectral regions in NGC 5044. 

As can be seen in Table [5] the iron abundance in EPIC is 
16% higher than the value deduced from RGS. The emission 
measure is lower by 24%, while the O/Fe ratio is 20% lower. 
Contrary to RGS, EPIC shows a significant temperature structure 
with the gdem model. Silicon, sulfur, argon and calcium to iron 
abundance ratios are high. The nickel to iron abundance is low. 

4. Discussion 

We accurately determined the metal abundances of the hot gas 
in two groups of galaxies. In the following sections we focus 
the discussion on four main points. First, the derivation of the 
temperature structure, including the fitting of the iron-L com- 



Parameter 


Kuo waem 


KLro gdem 


nrlC gdem 


Y (10 /u nT J ) 


5.45(17) 


5.45(16) 


4.44(4) 


W (keV) 


0.67(3) 


0.665(2) 


0.6713(13) 


cr T 


- 


< 0.012* 


0.057(4) 


l/a b 


<0.08 


- 


- 


P c 


0.65 


- 


- 


C/Fe 


1.7(4) 


1.7(4) 


- 


N/Fe 


1.8(3) 


1.8(3) 


- 


O/Fe 


0.72(5) 


0.73(4) 


0.61(3) 


Ne/Fe 


0.58(10) 


0.58(10) 


- 


Mg/Fe 


0.85(8) 


0.86(6) 




Si/Fe 






0.99(3) 


S/Fe 






1.09(5) 


Ar/Fe 






1.0(2) 


Ca/Fe 






1.9(5) 


Fe 


0.538(19) 


0.538(17) 


0.588(8) 


Ni/Fe 






0.54(9) 


C-stat / d.o.f. 


3864 / 2667 


3864 / 2667 


3601 / 1775 



Notes. Abundance ratios are given in proto-Solar units by Lodders 
(2003). Numbers between parentheses are the lcr errors on the same 
number of least significant digits. <0) Since l/a = 0, the wdem mean and 
peak temperatures are equal. <fc) 2a upper limit. (c) Kept fixed during the 
fit. 



plex and its effects on the reported abundance ratios. Second, 
we show how our abundance ratios compare to what was earlier 
reported in the literature. After this the supernova models and 
numbers using the abundances found in the gas are discussed. 
Finally we address the stellar origins of the carbon, nitrogen and 
oxygen abundances. 



4.1. Temperature structure 

In general, multi-temperature models fit the spectra of both 
groups well. In the EPIC spectrum of NGC 5 813 we find some 
system atic effects in the residuals near 1 .2 keV. Brickhouse et al.l 
(2000) encountered the same feature in Capella. They solve this 
problem by adding lines caused by transitions of Fe xvn, Fe xvm 
and Fe xix to the model. The abundances of Capella significantly 
increased when doing this. We recalculated the best fitted model 
using an experimental version of the line list that is being im- 
plemented in SPE^Q, which contains many more transitions, in- 
cluding the ones probably causing the feature. Doing this did 
improve the fit near 1.2 keV. To date it is however not possible 
yet to use the new list for fitting yet. 

In RGS, we find that some iron lines in the iron-L com- 
plex are not well fitted in both groups. In Table [6] we list 
the ten strongest lines between 15.3 A and 16.7 A in the fitted 
models to our sources. The eight brightest lines are the same 
for both sources but their order in the list is slightly different. 
Most of these strongest lines are due to the Fe xvm ion. We in- 
vestigated whether these model errors could be caused by the 
use of an outdated ionisation balance. We fitted t he spectra of 
NGC 5044 again using the ionisation balance by iBrvans et al.l 
(2009 j) instead of the de fault ionisati on balance of SPEX 
(lArnaud & Raymond! 1 19921 for iron and lArnaud & Rothenflugl 
1 19851 for the other elements). Statistically, the fit did not improve 
by this (C-stat 2824 cf. 2418 for the original balance). The fit 
around the 16 A line did not improve either. We also recalcu- 
lated the best fitted model using the forementioned experimental 
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Table 6. The ten strongest modeled lines around the 16 A blend, 
using the fitted model to the sources. 



Ion 


Wavelength 


Emission 




A 


xlO 47 photons s -1 


NGC 5044 


Fexix 


15,358 


11.7 


Fexix 


15.364 


16.1 


Fe xvm 


15.831 


24.3 


Fe xvm 


16.165 


26.1 


Fexix 


16.282 


26.4 


Fe xvm 


15.628 


38.9 


Fe xvm 


16.002 


41.0 


vm 


16.003 


62.1 


Fexix 


16.106 


87.2 


Fe xvm 


16.078 


120.9 


NGC 5813 


Fexix 


16.282 


7.7 


Fe xvm 


15,874 


7.9 


Fe xvn 


15.456 


9.9 


Fe xvm 


15.831 


16.9 


Fe xvm 


16.165 


17.3 


O vm 


16.003 


18.0 


Fexix 


16.106 


25.3 


Fe xvm 


15.628 


26.9 


Fe xvm 


16.002 


28.6 


Fe xvm 


16.078 


84.3 



5813 shows a single-temperature spectrum in RGS. This is fully 
consistent between the wde m and gdem fits of the RGS data 
and was previously found by Wer ner et al.l (12009b . In the cross- 
dispersion direction of RGS we do measure a significant but rel- 
atively small rise of the temperature near the outermost bin (be- 
tween 1' and 2')- In the EPIC data, the gdem models have a cr x 
of 0.057±0.004, corresponding to a 0.09 keV spread in the tem- 
perature. 

The analysis of Chandra data by Rand all et alJ (1201 ll) shows 
a high-resolution temperature map of the inner region of NGC 
5813 (r<3')- The temperatures range from ~0.6 keV to ~0.7 
keV. To investigate whether some, maybe unresolved, temper- 
ature structure would cause an effect comparable to our results, 
we simulated a spectrum consisting of two CIE components. We 
choose to set one of these components to have a temperature of 
0.6 keV and put the other at 0.7 keV. The normalisations were 
fixed to half the value found for NGC 5813 and all other param- 
eters were fixed to the best-fit values derived from the RGS data 
of NGC 5813. We then fitted this spectrum under the assumption 
of one single gdem distribution and found a crj of and a value 
for the ion to electron temperature ratio that is the same as we 
find in our fits. A fit to the NGC 5813 data using two distinct CIE 
models caused both temperatures to converge to the value found 
using a single gdem model. The effect found by leaving the ra- 
tio between the electron and ionisation balance temperature free 
could in both sources be caused by the part of the temperature 
structure that can not be resolved using the models. 



version of the SPEX line list. This also did not improve the fit 
near the blend. 

Work by lAudard et all (l200ll) using a MEKAL model on 
Capella, which has a thermal component with a temperature 
comparable to our sources, shows that the spectral lines around 
16 A are indeed not fitted well. The fact that the other line blend 
of Fe xvm and Fe xix is correctly fitted in one of the cases hints 
towards uncertainties in the unresolved spectral lines at these 
wavelengths. Since most of the lines of iron are fitted correctly 
when using the abundances that we find, the effect of single 
poorly fitted lines on the iron abundance should be limited. 
IWerner et al.l (2009) show an indication for resonant scattering 
of one of the Fexvn lines in NGC 5813 which could cause a 
slight underestimation of the iron abundances. Although they try 
to take this effect into account, the 16 A line is not perfectly fit- 
ted in that work as well. Since the blend at 16 A also consists 
of the O vm-Lya an resonant scattering effect on this line might 
also influence the underestimation of the flux at this blend. 

Line ratios that differ from the ones calculated using the ion- 
isation balance at a given temperature might be caused by non- 
equilibrium effects. To test for the presence of these effects, the 
CIE model allows for fitting of the ratio between the ionisation 
balance temperature and the electron temperature. The best-fit 
ratio we find is 1.21+0.03 for NGC 5044 and 1.30 ± 0.05 for 
NGC 5813. The fitting of this parameter improves the C-statistic 
with ~50. However, the problems with the 14 A and 16 A lines 
are not solved. Moreover, due to the high densities and colli- 
sional timescales in the cores of the groups of galaxies studied, 
it is highly unlikely that non-equilibrium effects are present. 

Clusters and groups of galaxies intrinsically contain mul- 
tiple gas components with different temperatures. Therefore, 
we need multi-temperature models to c onstrain temperature 
structures and abun dances correctly (e.g. iPeterson et al.l [2 003; 
iKaastraet alll2004l) . Knowing this, it is remarkable that NGC 



4.2. Abundances 

The oxygen abundances as measured by EPIC are lower than 
what is measured using RGS. Oxygen lines can not be well re- 
solved by EPIC due to its relatively poor spe ctral resolution and 
uncertain calibration (e.g. Ide Plaa et al.ll2007l) . Since we ignore 
EPIC-pn in the region of the oxygen line, instrumental effects are 
probably due to the cross-calibration between RGS and EPIC- 
MOS. Since the RGS spectra resolve the Ovm line at ~ 19 A, 
we use the RGS value. For NGC 5044 we could fit the EPIC 
spectrum leaving the instrumental power-law normalisation to 
vary within a factor 2. For NGC 5813, we needed to increase 
the factor to 4 and to leave the photon index break free. This 
effect could be due to a possible low contribution of quiescent 
soft-protons. We compared the abundances and temperature for 
those fits to a fit leaving only the normalisation free within a fac- 
tor of 2. The abundances of nickel, argon and sulfur did signifi- 
cantly differ between those fits. We therefore decided not to use 
the abundances of nickel, argon and sulfur for further analysis in 
NGC 5813. 

To check for a potential residual low-mass X-ray binary con- 
tribution to our spectra, we added a power-law to the spectral 
model in RGS. This did not improve the fit in any way, whether 
or not we fixed the photon index to a v alue of 1.5. For NGC 5813 
this is in agreement with the results of Rand all et al 

The values found in the literature for NGC 5044 and NGC 
5813 are summarised in Tables [7] and [8] Most other work 
has been carried out using spectra based on different extrac- 
tion regions, instruments or models. All these can have an ef- 
fect on the fitted temperatures and absolute abundance values. 
Abun dance ratios do howeve r generally not change with radius 
(e.g. [Simiones cu et al.l l2009). Because of this, the derived val- 
ues for those should be comparable. Differences in those values 
therefore are an indication of the systematic effects introduced 
by the usage of different instruments and models. 
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Table 7. Literature values of the parameters for NGC 5044 Table 9. Abundances reported in M87 bv lWerner et al.l (f2006). 



Parameter T03 B03 This work 

Instrument RGS ACIS" & EPIC RGS & EPIC 

r exlr Y*> 5' 2.5' b & 3' 

Plasma code 2 * CIE 2 * APEC gdem 

kT 1.07(3) 1.38(2) 0.910(5) 

0.804(4) 0.665(8) 

N/Fe 0.6(2) - 1.5(2) 

O/Fe 0.56(2) 0.42(4) 0.60(3) 

Ne/Fe 0.56(11) 0.75(19) 0.37(10) 

Mg/Fe 0.85(5) 0.87(3) 0.75(5) 

Si/Fe 0.97(7) 0.81(2) 0.788(12) 

S/Fe - 0.516(19) 0.843(16) 

Fe 0.74(3) 0.99(4) 0.616(16) 



Notes. T03: lTamura et al.U2003h. B03:lBuote et al.l d2003h . Abundance 
ratios in proto-Solar units by Loddersl (120031) . ( ^ on board of the 
Chandra satellite ( *' Half width of the region in cross-dispersion di- 
rection used for the analysis 



Table 8. Literature values of the parameters for NGC 5813. 



Parameter 


W09 


Rll 


Rll 


This work 


Instrument 


RGS 


RGS 


ACIS" 


RGS & EPIC 


'"extr 


0.25' b 


0.5'* 


2.9' 


2.5'* & 3' 


Plasma code 


CIE C 


2 * CIE d 


2 * APEC 


gdem 


kT 


0.645(8) 


>0.4 


0.671(11) 


0.665(2) 






<0.4 


0.35(3) 




N/Fe 


2.6(1.1) 






1.8(3) 


O/Fe 


0.71(15) 


0.80(14) 


0.27(6) 


0.73(4) 


Ne/Fe 


0.4(3) 


0.72(12) 


1.14(12) 


0.58(10) 


Mg/Fe 






0.83(7) 


0.86(6) 


Si/Fe 






0.95(7) 


0.94(2) 


S/Fe 






1.17(12) 


1.13(5) 


Fe 


0.75(9) 


0.81(8) 


0.71(4) 


0.538(17) 



Notes. W09: IWerner et al.l d2009b Rll: [R andall et aQ J201lh . 
Abundance ratios in proto-Solar units bv lLoddersl(l2003l) . ( '' ) on board of 
the Chandra satellite '*' Half width of the region in cross-dispersion di- 
rection used for the analysis <r) Model including four Gaussian compo- 
nents to correct for possibly present microturbulence effects. {d) Spectra 
were fitted using two CIE cooling flow models 



iTamura et all (120031) (T03) fitted the RGS spectrum of NGC 
5044 (obsID 0037950101) in the 0'-l' cross-dispersion re- 
gion using a mode l consisting of two temperature components. 
iBuote et all (l2003h (B03) used a two-temperature APEC model 
to constrain the parameters of NGC 5044 using Chandra. The 
abundance ratio of oxygen over iron is comparable between the 
two RGS observations and lower in Chandra. The nitrogen to 
iron abundance ratio we measure is higher than derived by T03 . 
This may be caused by the different background subtraction, ex- 
traction regions and temperature structure model (note that the 
background is relatively large at this wavelength and therefore 
sensitive to systematic effects). The same holds for the silicon 
value derived from RGS by T03. The value we derive from EPIC 
does however agree with the value of B03 derived from ACIS 
and EPIC. For the magnesium over iron ratio, the value we mea- 
sure is somewhat lower than found by both T03 and B03, but 
still acceptably close. The lower sulfur over iron ratio found by 
B03 can be explained by their higher temperature, at which the 
emission in the sulfur lines is higher and hence a lower abun- 
dances is needed to obtain the same line power. For neon, the 
values reported by the three studies do not agree well. Since the 
neon lines in the spectrum are relatively weak and blend with 



Parameter value 

kT m , im 1.82(4) = 

C/Fe 0.74(13) 

N/Fe 1.6(2) 

O/Fe 0.59(4) 

Ne/Fe 1.25(12) 

Mg/Fe 0.60(6) 

Fe 1.06(3) 

Notes. Abundance ratios in proto-Solar units by Lodders (2003). Based 
on fits to the data using a wdem model. 



iron, some systematic effects could still play a role because the 
often imperfect fit to the iron-L complex. We decided not to use 
it f or further analysis. 

IWerner et alJd2009b (W09) fitted the RGS spectrum of NGC 
5813 in a strip with width 0.5' (based on obsID 0302460101) 
in the cross-dispersion direction. For their fits, they take res- 
onant scattering of the strongest Fe xvn lines into account. 
iRandall etaf] (120111) (Rl 1) present results of an RGS fit to the 
same data set, using a strip of 1' and they also present the results 
of Chandra analysis. The temperature found by W09 is close to 
what we measure in the inner bins of the cross-dispersion pro- 
file. The oxygen over iron abundance ratio that we find is con- 
sistent with W09 and the RGS value found by Rll. The value 
reported by Rl 1 based on ACIS is much lower. As we noted be- 
fore, RGS has the resolving power to constrain the oxygen abun- 
dance using the O vm- Lya line as opposed to the CCD instru- 
ments where this line is part of a blend. The nitrogen over iron 
abundance ratio we find is consistent with the value by W09. 
There is a good agreement between the Chandra ACIS results 
and our work based on EPIC for the abundance ratios of magne- 
sium, silicon and sulfur over iron. However, the neon abundance 
is rather different and it is excluded from further analysis. 

We can compare the abundances that we find to those in other 
giant ellipticals in the cores of groups or clusters . The central 
galaxy of the Virgo cluster, M87, was studied by IWerner et al.l 
(2006). We show the values in Tableland compare our results 
with their findings because all three sources contain a giant el- 
liptical galaxy. The absolute iron abundance that we find is about 
half the iron abundance in M87. As mentioned before, absolute 
abundances are however sensitive to bias. The O/Fe ratio found 
for NGC 5044 is close to that of M87, while NGC 5813 has a 
somewhat higher value. The N/Fe ratio found for our groups is 
in good agreement with M87. We do find a value of C/Fe which 
is about a factor of 2 higher than in M87. It should however be 
noted that the model background of the RGS suffers from severe 
probl ems in the CCD ch ip on which the carbon line is projected. 

IWerner et al.l (1200 9) present the RGS-measured abundances 
of five bright elliptical galaxies, including NGC 5813. The fits 
are based on single-temperature models and include corrections 
for resonant scattering effects. The abundances are shown in 
Table [10] The fits presented in this table are based on relatively 
short observations (17 - 81 ks) which translates in relatively 
large error bars on the parameters. The ratio of nitrogen to iron 
is higher than solar for all objects. It appears that data from the 
shorter observations hint towards higher values for this ratio al- 
though we do not have enough data to confirm this. For NGC 
5813 the longer observation shows a lower abundance ratio. For 
the oxygen over iron ratio, the values of the abundances are all 
comparable with what we measure in our groups. As mentioned 
before, the neon abundance is very hard to constrain. The mea- 
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Table 10. Abundances reported by IWerner et al.l (2009). 



Parameter 


NGC 4636 


NGC 1404 


NGC 4649 


NGC 4472 


kT 


0.606(6) 


0.608(9) 


0.774(7) 


0.781(6) 


N/Fe 


2.5(6) 


3.4(1.3) 


1.5(9) 


1.6(6) 


O/Fe 


0.85(11) 


0.87(18) 


0.7(2) 


0.64(10) 


Ne/Fe 


0.60(16) 


1.2(4) 


1.5(5) 


1.4(3) 


Fe 


0.52(3) 


0.67(8) 


0.87(18) 


0.83(8) 



Notes. Abundance ratios in proto-Solar units by Lodders (2003). 



sured values are sometimes hardly significant and there is some 
spread. The absolute iron abundances do also show some spread. 
Those are however, as mentioned before, sensitive to biases. 

The abundances of carbon, nitrogen, oxygen (in RGS), sili- 
con and sulfur can be constrained based on single, isolated lines 
and are much more robust than the abundances of neon and mag- 
nesium because these metals are only constrained using lines 
that are blended with iron lines. The calcium and argon abun- 
dances measured are high with large error bars. Moreover they 
may be somewhat biased by residual background in the relevant 
energy band. The peak of the calcium and argon lines is 5 - 10% 
above the continuum, while the continuum shows systematic off- 
sets of 10 - 20% at the wavelengths of these elements (see e.g. 
Fig. [Tol l. In NGC 5044, the relative abundances measured with 
RGS are also consistent for the two DEM parametrisations. The 
iron abundance itself differs significantly, as well as the normal- 
isation of the hot gas model. 

The abundance ratios give an insight in the processes that 
create the metals in a cluster or group of galaxies. Most of the 
metals in galaxies are created by supernovae. Some metals, espe- 
cially carbon and nitrogen, are however not so efficiently created 
by supernovae and result from stellar winds of main sequence or 
AGB stars. We discuss both points in a bit more detail. 

4.3. Relative supernova yields 

For years, the oxygen and iron abundances of our Galaxy 
have been used to determine the ratio between type la and 
core-collapse supern ovae (see e.g. the reviews by lRanalll991t 
IWooslev et al.ll20(52h . In general type la supernovae are believed 
to yield most of the iron in the Galaxy while core-collapse su- 
pernovae would be responsible for the oxygen. The spread of the 
yields between the different models is however large. Attempts 
have been made to constrain the supernova ratios and models us- 
ing the abundances in the hot gas of clusters of galaxies (see e.g. 
Ide Plaa et alJl200l . 

The abundances that we measure in the gas do not directly 
reflect the chemical yield of supernovae during the total history 
of the group. But the abundances do provide constraints on the 
ratio of t ype la to core-collapse supernovae contribut ing to the 
gas (see lHumphrev & Buotd2006l:lde Plaa et alj2007h . If we as- 
sume that all iron and oxygen has been created by supernovae, 
we can use the ratio of these two abundances to find the ratio 
between type la and core-collapse supernovae that enriched the 
gas. We write 



sL-'(sL-<'-*(s) 



(5) 



SNcc 



Here X is the abundance of an element, Fe the iron abundance 
and / represents the fraction of the iron that was created by type 
la supernovae. The supernova models that we use are the best-fit 



models found bv lde Plaa et alJ (20071). For cor e-collapse super- 
novae, this is the model bv lNomoto et al.l (120061) . using s tars with 
Z = 0.02 and a Salpeter Initial Mass Function (IMF; ISalpeterl 
fl955l) . 

Since the yield of a core-collapse supernova depends on the 
mass, the typical core-collapse yield we use is the IMF weighted 
integrated yield (Tsuiimo to et al.l[T995l) . 

For the Type la supernovae, we use four different mod- 
els. Fjrs^_we_use two delayed detonation (DDT) models 
by iBadenes et al.l d2003l) . The DDTc model compared best 
to the abundances f ound in the Tycho supernova remnant 
(IBadenes et alj|2006l) . IBadenes et all d2008h pro pose that SNR 
0509- 67.5 is best fitted by the DDTa model from Ba denes et all 
(120031) . The difference between the DDT models is due to the 
different input parameters for the density of the ambient medium 
and the kinetic energy of the ejecta. The DDTa model predicts 
lower amounts of silicon to calcium than DDTc (Table [TO . 
We also compare our data to the W7 and WDD2 model from 
llwamoto et alJ (1 19991) because they are commonly used in the 
literature. These two models do also predict higher amounts of 
silicon and sulfur than the DDTa model. There is also some 
spread in the yields of argon and nickel between the four mod- 
els. The DDTc model predicts the highest argon value, which 
is still lower than we measure for NGC 5044. The nickel abun- 
dance predicted by the W7 model is higher than what we find 
while for the other models it is significantly lower. The amount 
of produced iron predicted by the DDTa model is higher than in 
the other models. This translates in the lower number ratio for 
comparable /. 

To find the value for /, we use the oxygen over iron abun- 
dance. We then use this value for / to predict the abundance ra- 
tios of the other elements expected by the models and compare 
these to our results. This is shown in Table [TTI We find a value 
fo r / which is consistent with the value found for NGC 5044 
by Bohrinaer et al. (2005), based on the abundances reported by 
Buote et al. (2003), using oxygen, iron and silicon abundances . 
The number is also consistent with M87 (IWerner et al.l 12006). 
This means that the the stellar population of our sources is sim- 
ilar to the population of M87. If we use the abundances of the 
solar neighbourhood as input (i.e. 0/Fe=l), we find a ratio of 
~20% of the supernovae being of type la. The stellar population 
of the group is significantly older than in the surroundings of 
our S un. By fitting super nova yields to multiple abundance val- 
ues, |de_PlaaetaL| (|2007) derived from a sample of clusters that 
44±5% of the contributing supernovae is of type la. This means 
that the stellar population in our groups seems to be slightly 
younger than or comparable to the average population in clusters 
of galaxies. The magnesium abundance predicted by the super- 
nova models is lower than the value we measured. The measured 
silicon abundances are higher than the values found using the 
DDTc model, and even more higher than the DDTa model. The 
DDTc model reproduces a high calcium abundance, but this is 
still lower than our measurements. 

We compare the abundances wi th the metallicities of stars 
in the disk of our own Galaxy (by IChen et al ] l2000t IShi etail 
120021: [Bensby & Feltzingl 120061) . Based on this we can indicate 
what the differences are between the large ellipticals in the cen- 
tre of groups and our Galactic environment. For comparison with 
Galactic stars, we use [X/Y] = logCY/F) - log(XYF) as pa- 
rameter, c onverting the abundances from the proto-Solar abun- 
da nces of Lodders (2003) to the photospheric solar abundances 
of lLoddersI (120031) . However, the difference between the pho- 
tospheric and solar abundance is 0.07 dex for all elements for 
which we present a ratio (C, N, O, Fe). The absolute values of 
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Table 11. Comparison of the Supernova yields to the measured 
abundances. 



Table 12. Comparison of the measured abundance ratios to the 
stars in the disk of our galaxy. 



Parameter 


DDTc 


DDTa 


WDD2 


W7 


measured 


NGC 5044 


/ 


0.84 


0.83 


0.83 


0.84 




C/Fe 


38% 
0.19 


31% 
0.19 


36% 
0.19 


39% 
0.21 


1.5 ±0.3 


N/Fe 


0.22 


0.23 


0.22 


0.22 


1.5 ±0.2 


O/Fe 


0.6 


0.6 


0.6 


0.6 


0.60 ± 0.03 


Mg/Fe 


0.47 


0.48 


0.49 


0.48 


0.73 ± 0.05 


Si/Fe 


0.67 


0.50 


0.74 


0.66 


0.788 ± 0.012 


S/Fe 


0.76 


0.51 


0.74 


0.62 


0.843 ±0.016 


Ar/Fe 


0.58 


0.36 


0.49 


0.38 


0.789 ± 0.013 


Ca/Fe 


0.96 


0.56 


0.71 


0.47 


1.86 ±0.13 


Ni/Fe 


1.01 


1.11 


1.16 


2.48 


1.60 ±0.05 



NGC 5813 



/ 


0.80 


0.79 


0.80 


0.80 






Ni a +Ncc 


32% 


25% 


32% 


33% 






C/Fe 


0.23 


0.23 


0.24 


0.25 


1.7 ± 


0.4 


N/Fe 


0.27 


0.27 


0.28 


0.27 


1.8 ± 


0.3 


O/Fe 


0.73 


0.73 


0.73 


0.73 


0.73 


±0.04 


Mg/Fe 


0.57 


0.57 


0.59 


0.59 


0.86 


±0.06 


Si/Fe 


0.73 


0.57 


0.81 


0.73 


0.99 


±0.03 


Ca/Fe 


0.98 


0.60 


0.74 


0.51 


1.9 ± 


0.5 



M87" 



/ 

Nig 
N, a +N„ 

C/Fe 

N/Fe 

O/Fe 

Ne/Fe 

Mg/Fe 



0.84 

38% 

0.19 

0.22 

0.6 

0.92 

0.47 



0.83 

31% 

0.19 

0.23 

0.6 

0.95 

0.48 



0.83 

36% 

0.19 

0.22 

0.6 

0.93 

0.49 



0.84 

39% 

0.21 

0.22 

0.6 

0.91 

0.48 



0.74 ±0.13 
1.62 ±0.21 
0.60 ± 0.03 
1.25 ±0.12 
0.60 ± 0.06 



Notes. Yields bas ed o n the models by 



Bade nes et all (2003), 



Iwam oto et al.l d!99S ) and iNomoto et al] d200€ ). Abundance ratios in 
proto-Solar units by Lodders (2003). The ratio of the iron created by 



supernovae (SN)Ia to the total amount of iron (see Eq. |5}> f, is cal- 
culated based on the O/Fe ratio using the assumption that all oxy- 
gen and iron has be en created by supernovae. <D) Values measured by 
Werner et al. (2006) are compared with the models used in this work. 



the abundances change, but the ratios are the same in both abun- 
dance systems. The abundance ratios measured in Galactic stars 
are calibrated on a measured Solar spectrum. We assume that this 
difference between the measured Solar spectrum and the recent 
Lodders abundances are small compared to our uncertainties. 

Since star formation has continued in our Galaxy, we can 
assume that the metallicities in young stars, which are located 
in the disk of the Galaxy, are comparable to the metallicities in 
the ambient medium. The oxygen to iron abundance ratio which 
we find for our sources is [O/Fe] = -0.222 + 0.022 at [Fe/H] = 
-0.210 + 0.011 for NGC 5044 and [O/Fe] = -0.137 + 0.018 at 
[Fe/H] = -0.269 + 0.014 f or NGC 5813. The ratios ar e low com- 
pared to the set of stars bv lBensbv & Feltzin g (2006). Therefore 
we can state that the contribution of core-collapse supernovae 
is relatively low in our sources compared to the core-collapse 
contribution in our Galaxy. 

4.4. Carbon and nitrogen progenitors 

For nitrogen and carbon, the measured abundance values are re- 
spectively factors 8 and 7 higher than predicted by supernovae 
only. For nitrogen, this num ber is much lower than the factor 
of 313 which was found bv lWerner et al.l (12006). The difference 
mostly rises from the fact that the supernova models we used for 



Ratio 


NGC 5044 


NGC 5813 


Typical range in disk stars 


[N/Fe] 
[N/O] 
[C/Fe] 


0.18±0.06 
0.40±0.06 
0.18±0.09 


0.26±0.07 
0.39±0.08 
0.23±0.10 


-0.2< [N/Fe] <0.2 fl 
-0.5< [N/O] <0.15" 
-0.2< [C/Fe] <0.2 fc 



Notes. Abundance ratios in solar units. (o) IChen et al] £2000) and 
Shietal] d2002l). < 6) Fro m both the thick and thin disk stars of 



iBens bv & Feltzing (2006) 



type la's include nitrogen creation. For the same reason the ra- 
tio we find for carb on is lower than the value of 10 reported by 
IWerner et al.l(l2006l) . 

Conclusions about the parent stellar population based on the 
carbon and nitrogen abundance in the hot intra-group medium 
should be drawn with care. First of all, these metals are still be- 
ing created and a fraction could still be locked up in stars. On the 
other hand, the stellar mass fraction of g roups and clusters is 10- 
20% of the gas mass (Ett oriet al.l l2009). This amount is compa- 
rable to the corresponding uncertainties on carbon and nitrogen 
and therefore the effect of those metals being locked up in stars 
is limited. Furthermore, the secondary creation of carbon and ni- 
trogen depends on the metallicity of the environment in which a 
star was formed. The formation of the less massive stars, which 
at this moment still enrich the gas, did however happen long 
ago during a relatively short period of the groups' life time. The 
metallicities of those stars around the time they formed are likely 
to be comparable to each other. In addition, our observations do 
not directly show how and when elements were ejected into the 
hot medium. As mentioned before, multiple mechanisms exist 
that enrich the medium with metals created in the galaxies. All 
mechanisms act on their typical time scale. Therefore, the mea- 
sured abundance ratios are the ratios in the intra-group gas, and 
not necessarily the ratios in the constituent galaxies. Since core- 
collapse supernovae, which can eject metals in a wind, seem to 
have decreased significantly in rate during the life time of the 
group, it is likely that other enrichment mechanisms took over. 
It was shown that jets from the AGN in the centre of the giant 
elliptica l galaxy can drag out metals from the galaxies into th e 
hot gas (ISimionescu et al]|2008l 120091; iKirkpatrick et al]|201 lb . 
Both NGC 5044 and NGC 5813 show evidence of AGN activity 
dDavid et aljl2009URandall et al.ll201 ll) . therefore AGN uplift is 
probably the currently dominating enrichment mechanism. 

The exact scena rio for the creation of nitrogen is no t yet fully 
understood (see e.g. lSchiavonll2007l:lPipino et al.ll2009l) . Models 
suggest that core-collapse supernovae produce much more oxy- 
gen than nitrogen. Nitrogen may instead be produced by low 
and intermediate mass stars. A high nitrogen over oxygen ra- 
tio would then hint towards lower present star formation. The 
comparison of the nitrogen to oxygen ratios of the two studied 
groups and the ranges of those values quoted in the liter ature 
(IChen et al.ll2000l IShi et al1l2002l: IBensbv & Feltzindl2006l) are 
shown in Table [T2l Based on the high nitrogen to oxygen abun- 
dance ratio in the two studied groups, we can conclude that fol- 
lowing the initial starburst the star-formation rate in these giant 
elliptical galaxies was significantly lower than in the Milky Way. 
This agrees well with the general assumption that the present star 
formation rate in elliptical galaxies is negligible. Due to the con- 
tinuous evolution of low and intermediate mass stars, nitrogen is 
still being enhanced while hardly any new oxygen is created. 
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When we compare the carbon and nitrogen to iron abundance 
ratios that we find to the typical abundance ratios in disk stars 
of our own Galaxy, we conclude that both the carbon and ni- 
trogen abundances are comparable to the stars with the highest 
ratios in disk of our Galaxy. The fact that the nitrogen and car- 
bon abundance ratios are higher than the ones in most of the 
stars of our Galaxy, while the abundance ratio of oxygen over 
iron is relatively low seems to point towards intermediate mass 
stars in AGB phase as sources of carbon and nitrogen creation in 
our sources, because the mass range of the progenitors of white 
dwarfs involved in type la supernovae (3-8 M Q ) is the same as 
for AGB stars. It should be noted that this conclusion is based 
on simple arguments and should be verified by more detailed 
population synthesis studies. 

5. Conclusions 

We present the results of the analysis of NGC 5044 and NGC 
5813 data, obtained with the RGS and EPIC instruments aboard 
XMM-Newton. From the analysis we conclude that: 

- The relatively low ratio of the oxygen over iron abundance, 
combined with the high abundance ratios to iron for carbon 
and nitrogen suggest that the bulk of the carbon and nitrogen 
are not created by massive stars which end their lives as core- 
collapse supernovae. 

- From the oxygen and iron abundance, and under the assump- 
tion of a Salpeter IMF, we find for both studied groups that 
about 80% of the iron is formed by type la supernovae. This 
translates to as roughly 30-40% of supernovae contributing 
being of type la while the other ~ 60-70% of supernovae 
then are core collapse. 

- The relative contribution of core collapse supernovae to type 
la supernovae is comparable with what was previously found 
for M87. This means that the giant elliptical galaxies in our 
study have a stellar and chemical evolution history compara- 
ble to the central giant elliptical galaxy of the virgo cluster. 
Both stellar populations tend to be older than in the surround- 
ings of our Sun. 

- A wdem and gdem model give comparable results in terms 
of abundances for the RGS observations of NGC 5044 and 
NGC 5813 

- Low temperature plasma models have problems fitting the 
line blends of Fe xvm, Fe xix and O vm near 14 A and 16 A. 
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